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Abstract 

Polymer-dispersed liquid crystal (PDLC) composites have long been a focus of study for their 
unique electro-optical properties which have resulted in various applications such as switchable 
(transparent/translucent) windows. These composites are manufactured using desirable “bottom- 
up” techniques, such as phase separation of a liquid crystal/polymer mixture, which enable pro¬ 
duction of PDLC films at very large scales. LC domains within PDLCs are typically spheroidal, 
as opposed to rectangular for an LCD panel, and thus exhibit substantially different behaviour in 
the presence of an external field. The fundamental difference between spheroidal and rectangular 
nematic domains is that the former results in the presence of nanoscale orientational defects in LC 
order while the latter does not. Progress in the development and optimization of PDLC electro- 
optical properties has progressed at a relatively slow pace due to this increased complexity. In this 
work, continuum simulations are performed in order to capture the complex formation and electric 
field-driven switching dynamics of approximations of PDLC domains. Using a simplified elliptic 
cylinder (microcapillary) geometry as an approximation of spheroidal PDLC domains, the effects 
of geometry (aspect ratio), surface anchoring, and external field strength are studied through the 
use of the Landau-de Gennes model of the nematic LC phase. 
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I. INTRODUCTION 


Traditional liquid crystal display (LCD) technology is principally based upon manipulat¬ 
ing optical properties of nematic liquid crystal (LC) thin films through the application of 
spatially localized electric fields. LCDs utilize the combination of a rectangular thin film ge¬ 
ometry, precisely controlled surface anchoring conditions, and the application of an electric 
field in order to manipulate a defect-free LC texture into a desired state DEI- A key aspect 
of traditional LCD technology is that conditions are engineered such that they do not impose 
topological constraints which result in the formation of orientational defects. Alternatively, 
an increasing number of different LC mechanisms and phases have been discovered which 
leverage the presence of orientational defects [3]. 

Polymer-dispersed liquid crystal (PDLC) films are one of the most studied LC-based ma¬ 
terials which involve LC dynamics with orientational defects present. One of the main appli¬ 
cations of these films is as switchable “privacy glass”, where an electro-optical property of the 
film is used such that an applied field drives the film to a transparent state and, upon release 
of the field, the film returns to a translucent state [1]. PDLC composites are formed through 
“bottom-up” manufacturing processes, mainly through photopolymerization-induced phase 
separation [T]. For low volume fractions of the LC component, the resulting composite 
morphology involves a spheroidally-confined dispersed LC guest phase in a polymer ma¬ 
trix phase. This spheroidal confinement imposes topological constraints on the dispersed 
LC domains which require the formation of orientational defects. Experimental observa¬ 
tion of these PDLC films in the presence of external fields has shown a rich and complex 
range of electro-optic behavior j3] depending on PDLC composition, chemistry, structure, 
polymer/LC anchoring conditions, and external field strength. 

Experimental studies have indicated that non-spherical droplet shape, specifically aniso- 
metric shape, is a key factor in PDLC film relaxation time following release of the external 
field | 3 ]. However, direct experimental observation of LC internal structure and dynam¬ 
ics for PDLC composites is challenging due to the length (nm — pm) and time (ns — ps) 
scales involved in LC dynamics. Simulation studies [5H22], however, have shed light on a 
far more rich landscape of internal structure than what is observable through experimenta¬ 
tion. Over the past two decades, simulation-based analysis has been used with increasing 
success, mainly focusing on cylindrical and spherical domains. Lattice-based simulations 
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[SHE] are able resolve sub-micron domains and have mainly been applied to study the effects 
of sub-micron cylindrical confinement where geometry and anchoring affects the stability of 
the nematic phase. However, continuum simulations [I2H21! have been able to overcome the 
length and timescales required to simultaneously capture defect dynamics (nanoscale) and 
domain shape (> pm). 

Continuum simulations of confined LC domains have been conducted using either 

Frank-Oseen vector theory [2S] or Landau-de Gennes tensor theory [MI]. While many of 
these past studies have focused on circular/spheroidal and elliptic/ellipsoidal [I2H20| 

confined LC domains, they have relied on Frank-Oseen theory which cannot capture orien¬ 
tational defects and phase transition. More recently, simulations of nematic LC confinement 
have been performed using the high-descriptive Landau-de Gennes tensor theory [211 122| , 
but these studies were limited to cylindrical domains and in the absence of an external field. 

In this context, the overall objective of this study is to predict both the formation and 
electric field-driven dynamics of nematic elliptic cylinder domains. While this geometry 
is a poor approximation of spheroidal and ellipsoidal domains observed in PDLCs 0, it 
has direct relevance to the study of nematic-filled capillaries which are of interest for fiber 
optics-based devices [26]. Furthermore, even though two-dimensional elliptic domains are 
not topologically equivalent to three-dimensional ellipsoidal domains, this “coarse” geometric 
simplification has been used in almost all past simulation-based work in the area except for 
in Ref. [20]. The Landau-de Gennes tensor model for the nematic phase is used in order to 
capture the presence of orientational defects, experimentally-relevant anchoring conditions, 
and phase transition. The specific objectives are to study the effects of geometry (aspect 
ratio) and anchoring conditions on: 

• the formation dynamics and equilibrium textures of nematic elliptic cylinder domains. 
This is similar to the state of the PDLC after quenching the phase-separated film. 

• the external field-driven dynamics of nematic elliptic cylinder domains. This similar 
to the state of the PDLG after application of an electric field to induce the transparent 
state. 

• the relaxation or “restoration” dynamics of nematic elliptic cylinder domains following 
release of the external field. This is similar to state of the PDLC after release of the 
electric field to return to the translucent state. 
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This study is performed with a few important assumptions; the hydrodynamic effects and 
thermal fluctuations of nematic order are neglected. Additionally, heterogeneous nucleation 
is assumed to be the dominant mechanism for formation of the domain from the disor¬ 
dered/isotropic phase, which is based on recent experimental observations [27j. 

The paper is organized as follows: first the theoretical and numerical bases for the simula¬ 
tion method are presented, next the nematic domain visualization and quantification meth¬ 
ods are described, then the results are presented and discussed for formation and switching 
dynamics, and conclusions made. 


II. BACKGROUND 

In this section the theoretical model and simulation methods are briefly described. The 
theoretical model is based upon the Landau-de Gennes theory for the nematic phase. This 
model is a Landau expansion for the nematic phase [2H] with respect to a tensor order 
parameter which enables simulation of phase transition and orientational defects (disclina- 
tions). Surface anchoring is modelled using a Landau-type expansion for the surface free 
energy [29l Chap. 4], as well. The governing equations are solved using the finite element 
method with second order implicit adaptive time-stepping EDI. Finally, visualization of 
the resulting tensor fields uses hyperstreamlines EH with a recently introduced topological 
seeding method [H2]- 


A. Landau—de Gennes Model 


The Landau-de Gennes model for the nematic free energy |2H] is expanded in terms of a 
second-order symmetric traceless tensor order parameter which is related to the anisotropic 
part of the magnetic susceptibility tensor |2H] • This tensor is known as the Q-tensor and is 
defined as: 


Q = XiTin + X2'mm + X3II 


S ( nn 



+ P {mm 


ll) 


( 1 ) 


where n/m/l are eigenvectors of Q which characterize the axes of molecular orientation, 
Aj are the eigenvalues of Q, S = |Ai is the uniaxial scalar nematic order parameter, and 
P = I (^2 — A 3 ) is the biaxial scalar nematic order parameter. S and P quantify the 
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extent to which the molecules conform to the local orientation where 5 ” = P = 0 in the 
isotropic phase and S' 7 ^ 0 , P = 0 in the uniaxial nematic phase. However, certain boundary 
conditions/geometries, external fields, and the presence of disclination defects can result in 
simultaneous order in multiple directions, resulting in biaxial order where both S and P are 
non-zero. 

The Landau-de Gennes model is able to capture both phase transition and elasticity 
through including terms in the expansion involving both Q and VQ: 


/,(Q, VQ, E, T) = T) + /„e(Q, VQ) + /^(Q, E) ( 2 ) 


where fb is the Helmholtz free energy density of the domain, fiso is the free energy contri¬ 
bution of the isotropic phase, f^b is the free energy contribution of bulk nematic ordering 
(thermodynamic), fne is the free energy contribution of the nematic distortions (elastic), /e 
is the free energy contribution from the electric field, and E is the electric field vector. The 
bulk nematic contribution is (Einstein notation): 

fnb 2^^^^ Pii)(QijQji) jk}Qki P ( 3 ) 

where aQfbfcsive material parameters and Tni is the theoretical second order isotropic/nematic 
transition temperature. The nematic elastic contribution is |HH] : 

fei = 2^1 (diQjkdiQkj) + 2^2 {d,Q^,dkQkj) + {dkQijd.Qik) ( 4 ) 


where L* are elastic material parameters which can be related to the fundamental modes of 
nematic deformation: splay, twist, and bend. In this study, the one-constant approximation 
is employed, making the L 2 and L 3 terms vanish. The electric field contribution is 

eii -|- 2e I 


fe = 
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where ey and e± are the dielectric constants parallel and perpendicular to the director n 
respectively. The above free energy formulation may be used to approximate the total free 
energy functional of a nematic domain in the presence of an electric field: 


F[Q] = [ hdV ( 6 ) 

Jv 

In order to simulate dynamics, with the previously mentioned assumptions of neglecting 
flow and thermal fluctuations, the dynamic equation used corresponds to the time-dependent 
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Ginzburg-Landau model or so-called Model A dynamics 
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where F = is the rotational viscosity of the nematic phase, and []'^^ is the symmetric- 

traceless component of the expression. At equilibrium this expression is equal to the Euler- 
Lagrange equation for the total free energy functional (eqn. |^. Expanding the functional 
derivative in the dynamic governing equation results in |33]: 
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B. Surface Anchoring Energy 


The two major types of physical anchoring conditions for nematic interfaces are homeotropic 
and planar, corresponding to anchoring parallel and orthogonal to the surface normal, re¬ 
spectively. In this work, homeotropic anchoring is utilized in that it uniquely constrains the 
preferred anchoring direction, whereas the use of planar anchoring would introduce addi¬ 
tional degrees of freedom in the boundary conditions. As with the bulk nematic domain, a 
surface free energy density can be formulated in terms of the Q-tensor field and the surface 
unit normal k |29j : 

fs fs,iso oikiQijkj (9) 


where fs,iso is the isotropic contribution to the surface free energy, a is the anchoring 
strength, and only terms up to first-order are retained. The total free energy (eqn. 
must now include a contribution from the bounding surface: 


F[Q]= [ hdV+ [ hdS 

Jv Js 


which imposes a boundary condition on the governing dynamic equation: 

= 0 


,, dfb 

+ kk- 


dQ 




'd{dkQij) 


( 10 ) 


( 11 ) 


C. Numerical Methods and Simulation Conditions 


Simulations were performed in two-dimensional elliptic cylinder geometries using the 
method of lines, where spatial discretization is achieved using the finite element method 
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TABLE I: Material parameters for 5CB. 


T 

m 

307.2 

K 

ao 

1.4 X 10*^ 

J/m^ K 

b 

1.8 X 10® 

J/m3 

c 

3.6 X 10® 

J/m3 

Li 

3.0 X 10“^2 

J/m 

e|| 

17 

(relative) 
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(relative) 


[HO] and time-stepping through a time-adaptive second-order implicit method. Mesh- 
independence simulations were first performed, finding that a uniformly distributed node 
density of 2.06 x lO^nodes/pm^ was sufficient using an error tolerance of 10“^°. Verifi¬ 
cation of simulation results was performed through convergence tests of equilibrium solu¬ 
tions. Simulation results were visualized using hyperstreamlines EH in conjunction with a 
topologically-informed seeding method for alignment tensor fields [H2]- Material parameters 
were used that approximate pentyl-cyanobiphenyl (5CB) [HSl ES] shown in Table 

Two types of simulations were performed: formation dynamics and switching dynamics. 
Formation dynamics simulations correspond to cooling of the PDLC domains from isotropic 
to nematic in the absence of an electric field; starting above the bulk nematic transition 
temperature (308.4 K) and cooling below to 307 K. Heterogeneous nucleation of the nematic 
phase was assumed to be the dominant nucleation mechanism based on recent experimental 
observations [2Z| . The initial conditions for these simulations assume a boundary layer that 
is uniaxial and well-aligned with the preferred orientational axis of the surface k{9), where 
0 is the polar angle of each surface point (see Appendix]^. Switching dynamics simulations 
correspond to application of an electric field to a fully-formed nematic domain. Simulations 
were performed with varying electric field strengths and with the field oriented along the 
major axis of the elliptic cylinder domain. The nondimensionalized form of the governing 
equation was used, and thus time is reported here as a dimensionless quantity, 

~ t Hr 

t = - , T = -— 

T ttQlni 
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III. RESULTS AND DISCUSSION 


Characteristic length scales were used to determine physically relevant domain sizes, 
surface anchoring energy, and electric field strengths for simulations. The most fundamental 
length scale is that of nematic ordering itself, which is determined through the competition 
between bulk thermodynamic and elastic contributions to the free energy [2H]: 




Li 


% (T - T„i) 


(13) 


which is the nematic coherence length. This length scale approximates the thickness of the 
nematic/isotropic interface. The nematic coherence length is typically ~ 4nm in the 
majority of LCs used in display technology. 

The nematic domain size was chosen to be on the order of 1pm, much larger than A„, 
with a constant area of 0.8 pm^ for all simulation geometries. For elliptic cylinder domains 
this constrains the cross-sectional area such that irab = 0.8 pm^ where a is the semi-major 
axis and b is the semi-minor axis of the ellipse. The pm-scale for LC domains is relevant to 
PDLCs used for privacy glass and other scattering applications [1]. Spheroidal geometries 
have been observed by Drzaic and others [Ij that are relatively complex, but range in aspect 
ratio from approximately R = (1,2]; thus this range of aspect ratios was used in simulations. 

Another characteristic length quantifies the competition between the electric field and 
the nematic elastic forces: 


A.= 



Li 


(14) 


- e^)E^ 

which is the dielectric coherence length [22] • The range of electric field strengths that was 
studied was [0Vpm“^,5 Vpm“^] which are typical field strengths for PDLC devices [2]. 
These field strengths result in Ae —>■ 10 nm which corresponds to a moderately strong electric 
field, but not to the extent that it could induce melting of the nematic phase. 

A similar characteristic length can be determined for the competition between surface 
anchoring and elastic forces: 


A.= 


h 

a 


(15) 


which is the surface extrapolation length [29]. As A^ —)■ 0 the surface anchoring effects 
dominate and the nematic alignment at the boundary governs the bulk texture; this is so- 
called “strong” anchoring. As A* —?• (X) the bulk nematic elasticity effects dominate and 






the nematic alignment at the surface is governed by the bulk texture. A value for the 
surface anchoring strength of 5 x 10“^ J/m^ was used with A* ~ 100 nm which corresponds 
to “weak” surface anchoring. 


A. Nematic Domain Formation and Equilibrium Texture 

Fig. El shows visualizations of the Q-tensor field for nematic elliptic cylinder domains 
with aspect ratios R = (1.05,1.6,2}. The dynamics observed in these simulations are rep¬ 
resentative of those for all aspect ratios studied, i? = (1,2]. A sequence of three distinct 
growth regimes was observed: free growth, interface impingement/defect formation, and 
bulk relaxation. The free growth regime involves the stable nematic phase growing into the 
unstable isotropic phase such that the texture is approximately commensurate with the an¬ 
choring conditions. The interface impingement/defect formation regime follows, where the 
nematic/isotropic interface impinges on itself resulting in the simultaneous formation of a 
pair of orientational defects along the major axis. For all simulations, the type of orienta¬ 
tional defects, or disclinations, observed were wedge-type with strength [371 Chap. 2]. 
These disclinations are formed in order to resolve the topological constraints imposed by the 
confinement geometry and anchoring conditions. Finally, the bulk relaxation regime follows 
impingement where the fully formed nematic texture relaxes to its equilibrium state through 
simultaneous disclination motion towards the ellipse focal regions and bulk reorientation. 

The growth and formation mechanism for the i? ~ 1 case has been studied in past work 
for two-dimensional domains Uni ED El ES]. Rey and Sharma developed a texture phase 
diagram using the Landau-de Gennes model for the circular geometry [22]. They predict 
the splitting of a -|-1 disclination into a pair of -|-^ disclinations located at the center of 
the domain for domain sizes on the order of A„. Given that the simulated domains are 
micron-scale, the defect splitting observed in the R ~ 1 case is in agreement with these 
past results. The resulting domain texture is uniformly oriented in a region parallel to an 
axis and with the pair of disclinations located along the axis orthogonal to it. Additionally, 
simulation results are found to also be consistent with past work [25] predicting a transition 
from textures with disclinations to an “escape” texture, absent of defects, under certain 
conditions (anchoring strength, domain size, etc). The growth and formation mechanism 
for the elliptic case i? > 1 is observed to be significantly different and does not exhibit 
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FIG. 1: (Color online) Hyperstreamline visualizations of the Q-tensor fields during the 
formation process for nematic domains: (a) R = 1.05 with t = 50, 1.82 x 10^, 2.43 x 10^, 
4.63 X 10^ 6.40 X 10®; (b) R = 1.6 with t = 50, 1.35 x 10^ 1.93 x 10^ 4.02 x 10®, 
4.56 X 10®; and {c) R = 2 with t = 50, 1.19 x 10^ 1.71 x 10^ 3.09 x 10®, 3.35 x 10®. 
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disclination splitting, even though the topologies of the elliptic and circular domains are 
equivalent. 

The underlying formation mechanism for elliptic domains is well-explained by the approx¬ 
imation of Wincure and Rey j3H] of the velocity of a uniaxial isotropic/nematic interface. 
They have shown that its velocity v is proportional to the difference in energy between the 
nematic and isotropic phases AF and capillary force C [38] : 

/3v = C-AF (16) 

where /5 is an effective viscosity term. During the free growth regime, in every case (Fig. [^, 
as the isotropic/nematic interface approaches the center of the domain, the capillary force 
grows inversely proportional to the radius of the central isotropic region. The first deviation 
of the elliptic formation mechanism is observed due to this competition of forces. For a 
circular domain, v{6) is essentially constant but for the elliptic domain it varies from a 
maximum at the interface regions closest to the minor axis to a minimum at those regions 
closest to the major axis. 

As the radius of curvature of the interfaces closest to the major axis approaches a critical 
value proportional to the nematic coherence length A„, the capillary force approaches the 
difference in free energy driving force {C —>■ AF). This results in a critical slowing down 
of the interface and the transition from the free growth regime to the impingement/defect 
formation regime. The conditions under which this occurs for the elliptic case cannot result 
in the formation of a -|-1 disclination which later splits, as is observed for circular domains 
[22] . Instead, a pair of -|-^ disclinations form directly along the major axis of the elliptic 
domain near the two (ellipse) focal points. Simultaneously, the isotropic/nematic interfaces 
in the central region impinge, forming a well-aligned central region along the minor-axis of 
the elliptic domain. This defect formation mechanism is the confinement-driven analogue 
to the defect “shedding” mechanism discovered by Wincure and Rey for free growth of 
nematic droplets in an isotropic matrix phase [36]. In both cases, defects are formed at the 
isotropic/nematic interface due to frustration between bulk droplet texture and interfacial 
anchoring. 

Finally, the fully nematic domain relaxes towards the equilibrium state, shown in the 
final sets of images in Fig. for the domains with aspect ratios R = {1.05,1.6,2}. This 
relaxation involves simultaneous motion of the disclinations along the major axis and bulk 
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reorientation. The relaxation mechanism for circular domains has been shown to be governed 
by the competition of bulk elasticity and surface anchoring strength [22] • Elliptic domains 
introduce an additional contribution: the variation of the curvature of the boundary. Fig. 

shows the evolution of the distance between defect cores for the formation process of 
each of the domains simulated. At equilibrium, the defect separation distance is found to 
increase with increasing aspect ratio which occurs without any changes in bulk elasticity 
or surface anchoring strength. This behavior can be explained by quantifying the mean 
curvature imposed by the elliptic boundary conditions. Figs. |3aH3b| show the schematic of 
an ellipse in polar coordinates and the mean curvature k as a function of 9 for ellipses of 
increasing aspect ratio. As aspect ratio increases, the curvature of the boundary regions 
increasingly becomes distinct: (i) a high-curvature region outward by the ellipse focal points 
and (ii) low-curvature regions elsewhere. Thus as aspect ratio increases the combination 
of geometry and anchoring effects increasingly impose highly localized deformation of the 
nematic in the focal regions of the ellipse, which results in equilibrium textures with defects 
located in these regions. 



FIG. 2: (Color online) Evolution of the distance d between defects for different aspect ratio 
equilibrium elliptic domains. In order from bottom curve to top, the aspect ratio is 

i? = {1.05,1.2,1.4,1.6,1.8,2.0}. 


B. Field-switching Dynamics 


Field-switching dynamics simulations were performed for the range of aspect ratios R = 


(1,2] using the results from Sec. Ill A as initial conditions. Electric field orientation was 
chosen to be parallel to the major axis which results in the most significant field-driven effect 
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0 (rad) 


(a) (b) 

FIG. 3: (a) Schematic of an ellipse in polar coordinates; (b) (color online) plot of the mean 
curvatnre k of the ellipse bonndary versus polar angle 9. 

on domain texture. The results of these simulations are presented and discussed by hrst 
focusing on the equilibrium domain textures during application of the electric field and after 
its release. Following this, the dynamics of the transition from the initial equilibrium domain 
to the field-driven domain and after release of the field are presented and discussed. In 
order to interpret the complex textures and texture dynamics resulting from these transient 
simulations, the droplet order parameter is introduced j39] , 

Qd = V-^ f QdV (17) 

Jv 

which can be further decomposed into the droplet director rid and droplet scalar order pa¬ 
rameter Sd- In all elliptic cases the droplet director rid is initially parallel to the minor axis 
due to the homeotropic anchoring conditions. The magnitude of Sd relates to the optical 
properties such that as Sd decreases, the domain more efficiently scatters light and as it 
increases, the domain more efficiently transmits light. 

1. Field-driven Equilibrium Textures 

Simulations of application and release of an electric field, with strengths ranging from 
(0Vpm“^, 5Vpm“^] were performed. Two different field-switching regimes were observed 
where the droplet director rid either remained constant (low field strength) or reoriented 
(high field strength). Visualizations of the equilibrium Q-tensor fields for both regimes are 
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shown in Fig. and Fig. respectively. These results support the assumption that a 
critical field strength Ec exists depending on the domain aspect ratio, anchoring strength 
and LC material properties. 

For the E < regime in the field-driven state, there is a small response in the domain 
texture resulting in a change in the droplet scalar order parameter Sd, but not the droplet 
director. Defects are driven inwards along the major axis which results in larger field-aligned 
focal regions; this mechanism similar to the domain texture prior to the bulk relaxation 
regime observed during formation. Upon release of the field the domain is restored to the 
initial equilibrium texture resulting from the formation process. 

For the E > Ef. regime in the field-driven state, there is a large response of the domain 
where disclinations transition to be aligned along the minor axis and the domain becomes 
strongly field-aligned through reorientation of n^. This is achieved through the domination 
of electric field forces over surface anchoring forces, although defect “escape” is not observed 
as it could be for very strong electric field or very weak surface anchoring conditions [30]. 
Upon release of the field the domain is restored to the initial equilibrium texture resulting 
from the formation process. 



0.17 0.20 0.23 0.27 0.33 

FIG. 4: (Color online) Hyperstreamline visualizations of field-driven textures with E < Ec 
for domains with aspect ratio (a) R = 1.05, (b) R = 1.6, and (c) i? = 2.0. 


Fig. I shows droplet scalar order parameters Sd for simulated equilibrium domains at 
different field strengths both below and above Ec- As electric field strength is increased (but 
still held below Ec) the effect is minimal on Sd- This can be explained through observation 
of the field-driven domain textures in Fig. while the area aligned with the field (focal 
regions) increases in size and uniformity, this is achieved through a simultaneous reduction 
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0.17 0.20 0.23 0.27 0.33 

FIG. 5: (Color online) Hyperstreamline visnalizations of field-driven textnres with E > E^. 
for domains with aspect ratio (a) R = 1.05, (b) R= 1.6, and (c) R = 2.0. 


in the central region that is aligned orthogonal to the field. Thus the net increase in nematic 
order in the focal regions is almost completely offset by decreased order in the central region. 

As electric field strength is increased above Ec there is a significant deformation of the 
domain texture, both with respect to Sd and n^. This response is similar to the Fredericks 
transition observed in planar LC domains exposed to an external field [2H] • In all simulations 
the droplet director rid is observed to reorient parallel to the applied field vector (major axis). 
The observed response of Sd after application of the field is more complex. For low aspect 
ratios i? < 1.5 the electric field is found to result in a field-aligned nematic domain with 
droplet scalar order parameter greater than in the absence of the field. This is the traditional 
mechanism associated with PDLC films used for privacy glass, where in the “off” state they 
scatter substantially more than in the “on” state. For higher aspect ratios i? > 1.5 the 
opposite behavior is observed; even though the domain texture is field-aligned, the droplet 
scalar order parameter is substantially lower in the field-driven state compared to without 
the field. This corresponds to an increase in light scattering, regardless of the orientation 
of the droplet director rid- This “reverse-mode”-like mechanism is typically associated with 
PDLC films formed using a nematic LC with negative dielectric anisotropy m- 


2. Field-driven Dynamies Textures 

Fig. 0 shows simulation results of the field-driven dynamics for the R = 2 case for E = 
5V pm“^ > Ec- These results are representative of simulation results for all domains R = 
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FIG. 6: (Color online) Plot of droplet scalar order parameter Sa versus electric field 
strength E for R = (1,2] where o and • glyphs correspond to _L and || E, 

respectively. 


(1,2]. A sequence of three distinct dynamic regimes were observed: bulk growth/recession, 
disclination/bulk rotation, and bulk relaxation. 

The bulk growth/recession regime involves the simultaneous growth of the field-aligned 
focal regions and recession of the field-orthogonal central region which involves the motion 
of the disclinations inwards along the major axis. This regime was observed for both the 
E < Ec and E > E^, cases, where in the former the net change in the domain texture was 
minimal (see Fig. [^. In the latter case (E > Ec) as shown in Fig. 9c, this regime involves a 
monotonic decrease in the droplet scalar order parameter —>■ 0, corresponding to a radial 
texture (maximal light scattering). 

The disclination/bulk rotation regime follows the bulk growth/recession regime; as the 
distance between disclination defects decreases (along the major axis) the repulsive nematic 
elastic forces approach that of the applied field. At that point, the defect separation distance 
becomes constant and rotation occurs. The initiation of this regime is found to involve a 
sharp transition of the droplet director from being orthogonal to parallel to the field direction 


(Fig. 9c). This discontinuous transition of is enabled by the radial texture of the droplet, 
where Sd — Q making a degenerate quantity. Disclinations and bulk nematic texture 
simultaneously rotate about the center of the domain, increasing domain field alignment 
rapidly (Fig. [^ . 

The bulk relaxation regime involves simultaneous rotation and expulsion of the disclina- 
tions from the central region along the minor axis. The rotation process ceases as disclina- 
tions approach the boundaries and surface anchoring forces balance out bulk elastic forces. 
Higher electric field strengths likely exist which would overcome surface anchoring and re- 
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suit in an “escape” texture m, but they exceed the field strengths typically used in PDLC 
applications [2], the focus of this study 

The observed mechanism predicted by simulations provides a more refined understanding 
of the mechanism predicted by Drzaic [Ij. As opposed to short timescale bulk reorienta¬ 
tion followed by long timescale defect motion, simulations predict that the mechanism is 
instead through the simultaneous motion of disclinations, growth of field-aligned regions, 
and recession of nonaligned regions. As originally predicted by Drzaic, the long timescale 
component of the mechanism is through the motion of defects transitioning from one axis 
of the domain to the other, but the short timescale mechanism also involves linear defect 
motion and reorientation dynamics in the outer regions of the domain. 



0.17 0.20 0.23 0.27 0.33 

FIG. 7: (Color online) Hyperstreamline visualizations of the field-on texture dynamics for 
the R = 2 domain with E = 5 V/pm and t = {0, 5.25 x 10^, 1.05 x 10®, 1.55 x 10®, 

1.89 X 10®, 2.11 X 10®, 2.33 x 10®, 2.71 x 10®, > 5.00 x 10®}. 

Fig. 1^ shows simulation results of the domain dynamics upon release of the field. The 
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regimes observed following release of the electric field are found to be similar to that of 
driven mechanism, but occurring in reverse. There are two significant differences: timescale 
and droplet director evolution. The timescales associated with each regime are an order of 
magnitude larger compared to their electric field-driven analogues, which is expected due to 
the different in magnitudes of the surface anchoring and electric field strengths. The surface 
anchoring strength governs the magnitude of the restoring forces which, in turn, govern the 
timescale for nematic orientation dynamics. They can be compared through their charac¬ 
teristic lengths, where ^ ~ 10 for the simulation conditions used. Additionally, during 
the disclination/bulk rotation regime the droplet director continuously rotates instead of 
exhibiting a sharp transition. This is due to the fact that the droplet scalar order parameter 
Sd does not decrease to zero during growth/recession regime. Thus a continuous rotation is 
required in order for the droplet director to return to the equilibrium configuration, orthog¬ 
onal to the field direction. Equilibrium textures in all restoration simulations were found to 
be the same as those from the formation simulations, indicating that under the simulation 
conditions the electric-field induced deformations of the nematic domain were viscoelastic. 



0.17 0.20 0.23 0.27 0.33 


FIG. 8: (Color online) Hyperstreamline visualizations of the release texture dynamics for 
the R = 2 domain, continued from Fig. (t = 0), with t = {1.79 x 10®, 6.75 x 10®, 
8.55 X 10®, 1.03 X 10^ 1.26 x 10^ > 1.60 x lOT 
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FIG. 9: (Color online) Droplet scalar order parameter evolution versus time for domains 
with a range of aspect ratios R = (1,2] with electric field applied (a,c) and following 
release (b,d): (a) E — 1 V/pm < Ec (driven); (b) E = 1 V/pm < Ec (release); (c) 

E = bY /pm > Ec (driven); and (d) E = bY /pm > Ec (release). Droplet director 
orientation is indicated by line type: solid and dashed correspond io na E E and 

Ud II E, respectively. 


IV. CONCLUSIONS 

In this work, a simulation-based study was performed on the formation and electric field 
switching dynamics of elliptic cylinder nematic domains. The observed nematic reorientation 
dynamics were found to have a complex dependence on geometry (aspect ratio), surface 
anchoring strength, and external field strength. Both formation and reorientation dynamics 
were found to be governed by the presence and motion of nematic disclination defects within 
the domain. Geometry of the domain, specifically aspect ratio, was found to strongly effect 
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domain texture by providing regions of high curvature to which nematic defects are attracted. 
Simulations also predict the presence of a geometry-controlled transition from nematic order 
enhanced by an external field (low aspect ratio) to nematic order frustrated by an external 
field (high aspect ratio). 

Equilibrium and dynamic behavior of elliptic nematic domains are found to significantly 
differ from circular ones, which opens up new possibilities for electro-optical mechanisms for 
applications of nematic-filled capillaries. Experimental validation of the presented results 
is needed, which can be achieved through comparison with future electro-optical dynamics 
measurements of nematic-filled capillaries as has been done for PDLC films jl]. Einally, the 
results presented support the use of simulation-based methods for rational design of PDLC 
optical functional materials tailored to application-specific requirements for nematic texture 
and switching dynamics. 
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Appendix A: Initial Conditions for Heterogeneous Nucleation 


The initial conditions for all simulations assume a boundary layer that is uniaxial and 
aligned with the surface normal k(9), where (r, d) are the polar coordinates of each surface 
point. A linear decay was used for the uniaxial nematic order parameter S (eqn. [^: 


Sinit{r,0) = Sb 


^ - rsjO) + 


for rs{9) - Xn < r < rs{9) 


(Al) 


S = 0 for 0 < r < rs{9) — 


(A2) 


where Sb is the bulk order parameter value at the simulation temperature, is the nematic 
coherence length (eqn. [^, and rs{9) is the radial coordinate of the ellipse surface. The 
expression for rs{9) is: 


c(^) 


a 

\/h‘^ cos^ 9 -\- a^ sin^ 9 


(A3) 
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where a/b is the length of the major/minor ellipse axis. The nematic director field (eqn. 
is assumed to be well-aligned with the ellipse surface normal: 

n = acos(0)ei + fesin(0)e2 (A4) 


resulting in the initial Q-tensor field: 


Q init 2 ^ 


(A5) 
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